% RigiditeElementaire retourne la matrice de rigidité elementaire associee
% a l'element element
% element : vecteur contenant les abcisses des noeuds de l'element
% poutre : structure contenant les donnees poutre. Il faut au moins un
% champs contenant le module d'Young (poutre.E) et un champ contenant le
% moment d'inertie de la section droite (poutre.I)

function Kel = RigiditeElementaire(element,poutre)

l = element(2) - element(1) ;

Kel = [6,3*l,-6,3*l;
       3*l,2*l^2,-3*l,l^2;
       -6,-3*l,6,-3*l;
       3*l,l^2,-3*l,2*l^2] ;

Kel = Kel * 2*poutre.E*poutre.I/(l^3) ;

end
